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Abstract. We examine the behaviour of an axisymmetric 
galactic dynamo model with a radial accretion flow in the 
disc. We also introduce a vertical magnetic flux through 
the galactic midplane, to simulate the presence of a large 
scale magnetic field trapped by the galaxy when form- 
ing. The trapped vertical flux is conserved and advected 
towards the disc centre by the radial flow. We confirm 
that accretion flows of magnitude several kms -1 through 
a significant part of the galactic disc can markedly inhibit 
dynamo action. Moreover, advection of the vertical flux 
in general results in mixed parity galactic fields. However, 
the effect is nonlinear and non-additive - global magnetic 
field energies are usually significantly smaller that the sum 
of purely dynamo generated and purely advected field en- 
ergies. For large inflow speeds, a form of 'semi-dynamo' 
action may occur. 

We apply our results to the accumulation and redis- 
tribution, by a radial inflow, of a vertical magnetic flux 
captured by the Galactic disc. Taking representative val- 
ues, it appears difficult to obtain mean vertical fields near 
the centre of the Milky Way that are much in excess of 
10 /iG, largely because the galactic dynamo and turbulent 
magnetic diffusion modify the external magnetic field be- 
fore it can reach the disc centre. 

Key words: Magnetic fields - MHD - ISM: magnetic 
fields - Galaxy: centre - Galaxies: ISM - Galaxies: mag- 
netic fields - Galaxies: spiral 



1. Introduction 

Dynamo systems only rarely occur isolated magnetically 
from their environment. To assume such isolation may be 
a plausible assumption for, e.g., stellar dynamos, but the 
situation is not so clearcut for spiral galaxies. Neverthe- 
less, it is generally supposed that the dynamo acts inde- 
pendently of any trapped vertical flux, and the interac- 
tion of a dynamo system with its magnetic environment 



has rarely been addressed. However, the effects of external 
magnetic fields have been invoked to explain the origins 
of various features of astrophysical discs, such as vertical 
dust filaments in spiral galaxies (Sofue 1987), strong verti- 
cal magnetic fields in the centre of the Milky Way (Sofue & 
Fujimoto 1987, Chandran, Cowley & Morris 2000) and, in 
a different astrophysical (but not conceptual) context, the 
launching and collimation of outflows in accretion discs 
around young stars and compact objects (Blandford & 
Payne 1982, Mestel 1999). It is widely accepted, however, 
that both accretion discs and the discs of spiral galaxies 
are sites of dynamo action (Parker 1979; Zeldovich, Ruz- 
maikin & Sokoloff 1983; Ruzmaikin, Shukurov & Sokoloff 
1988; Beck et al. 1996; Kulsrud 1999; Mestel 1999). In this 
paper we discuss, in a galactic context, the nonlinear in- 
teraction between a disc dynamo and a trapped vertical 
flux redistributed by the radial inflow of matter. 

The most important aspect of the problem in the galac- 
tic context is the origin of the vertical magnetic field in the 
Galactic centre where it can be as strong as 1 mG within 
the inner 200 pc radius, implying vertical magnetic flux of 
about 40 G pc 2 if the field is unidirectional and pervasive 
(Morris & Serabyn 1996). The field geometry is apparently 
dipole-like, in contrast to the quadrupole symmetry dom- 
inant near the Sun. With a gas density of 10-100 cm -3 
between molecular clouds and the turbulent velocities of 
20-30 kms -1 (e.g., Morris & Serabyn 1996), the field is 
100-10 times above the equipartition level with the turbu- 
lent kinetic energy. These features suggest that the vertical 
magnetic field in the Galactic centre might arise from an 
accumulation of an external magnetic flux advected from 
large radii over the Galactic lifetime (Sofue & Fujimoto 
1987, Chandran et al. 2000). However, the Galactic disc is 
magnetically active, being a site of dynamo action, and so 
it is not clear whether or not any vertical magnetic field 
can be dragged through it. Our results suggest that it is 
rather implausible that the vertical magnetic flux at the 
Galactic centre results from the redistribution of magnetic 
field trapped at an early stage of Galactic evolution. 
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2. Formulation of the problem 

The interaction of a nonlinear dynamo system with an 
external magnetic field can be non-trivial. Drobyshcvski 
(1977) proposed the idea of the 'semi-dynamo', a system 
which amplifies magnetic field at the expense of kinetic 
energy like a self-excited dynamo, but in which the field 
would decay without the presence of an external magnetic 
field. Convective dynamos in the presence of an imposed 
magnetic field have been studied by Sarson et al. (1997) 
and Sarson, Jones & Zhang (1999) in application to mag- 
netic fields of Io and Ganymede, Jovian satellites orbiting 
within the giant planet's magnctosphcre. These authors 
report on magnetic states with self-maintained magnetic 
field which, however, revert to a state of nonmagnetic con- 
vection if the external field is withdrawn or weakened. In 
this model, a uniform magnetic field parallel to the ro- 
tation axis is included as a source term in the induction 
equation, although the results remain the same when the 
external magnetic field is introduced through boundary 
conditions (G. R. Sarson, private communication). There 
is a plausible connection between the behaviour of these 
systems and that of Drobyshevski's semi-dynamos. 

There have been a number of studies of stellar dynamos 
with a fixed magnetic field as boundary condition at the 
bottom of the convection zone, simulating the presence 
of an interior field frozen in to the radiative core (e.g. 
Pudovkin & Benevolenskaya 1984, Boyer & Levy 1984, 
Moss 1996). The effect here was less fundamental, being to 
change properties such as parity, rather than cause semi- 
dynamo action. 

A two-dimensional mean-field dynamo with external 
magnetic field has been considered by Reyes-Ruiz & 
Stcpinski (1997) who, however, neglected any radial veloc- 
ity and obtained a solution of the nonlinear (a-quenched) 
dynamo equations with external magnetic field by merely 
adding a constant vertical field to the solution of the 
standard dynamo equations. The physical importance of 
their results for realistic degrees of nonlinearity is there- 
fore questionable. 

The inward advection (dragging) of an external mag- 
netic field by accretion flow in a Keplerian disc has been 
considered by Lubow, Papaloizou & Pringle (1994) (nu- 
merical solutions have been presented by Reyes-Ruiz & 
Stcpinski 1996). These authors disregarded any dynamo 
action and considered only the balance between advection 
of a purely poloidal magnetic field and radial diffusion. 
Their results show that efficient advection of an exter- 
nal magnetic field into a viscous disc is only possible if 
the magnetic diffusivity is much smaller than the kinetic. 
Otherwise inward advection is balanced by radial mag- 
netic diffusion and the steady-state magnetic field in the 
disc can only be weak. 

There are several possible types of external magnetic 
field for galaxies and accretion discs. This can be a field of 
cosmological origin. Current models of magnetic field gen- 



eration in the early Universe predict that the strength of a 
primeval field at the galactic scales can only be negligible 
(e.g. Beck et al. 1996), so this option has lost its impor- 
tance. For galaxies belonging to a galaxy cluster (but not 
the Milky Way), the intracluster magnetic field (Kronberg 
1994) can be a significant source of external magnetic flux. 
The strength of this field is about 1 /iG, but only a frac- 
tion of it may have the appropriate poloidal structure in 
the galactic frame. Away from rich galaxy clusters, the 
observational upper limit on the cosmological magnetic 
field (< 3 x 1(T 12 G, Sciama 1994, Beck et al. 1996) im- 
plies an upper limit of 3 x 10~ 10 G for a quasi-uniform 
magnetic field captured by and compressed in a collaps- 
ing protogalaxy (for radial compression by a factor 10). 
For accretion discs in an active galactic nucleus or in a 
stellar system, the interstellar magnetic field of the par- 
ent galaxy (~ 1 [iG for the total field of which a fraction 
would be poloidal field as seen by the disc) can serve as 
an external field. Stronger external magnetic fields origi- 
nating in the central star can occur in accretion discs of 
binary stellar systems (e.g. Bardou & Heyvaerts 1996). 

Although the relative orientation of the disc plane and 
the external magnetic field can be arbitrary, we neglect 
any magnetic field component parallel to the disc plane 
and allow only for a component parallel to the rotation 
axis. Any horizontal component will be twisted by dif- 
ferential rotation until the separation of magnetic line 
segments of alternating sign becomes small enough for 
the field to be destroyed by turbulent magnetic diffusion 
and/or reconnection (e.g. Sects. 3.8 and 3.9 of Moffatt 
1978). The only way to avoid or postpone the expulsion 
is by reducing the magnetic diffusion as, e.g., in Howard 
& Kulsrud (1997) and Chandran et al. (2000), who only 
consider the effects of a relatively weak ambipolar diffusiv- 
ity. Observations directly indicate the presence of intense 
turbulent motions in the interstellar gas of spiral galaxies 
(Minter & Spangler 1996). Even though interstellar mag- 
netic fields are dynamically important, there are no ap- 
parent indications of any modification of the turbulence 
that could result in any significant suppression of tur- 
bulent diffusion. Furthermore, Howard & Kulsrud (1997) 
neglect magnetic buoyancy, which can be efficient in re- 
moving horizontal magnetic field from the disc on a time 
scale of order the Alfven crossing time over the disc scale 
height. Magnetic buoyancy possibly presents insurmount- 
able difficulties for a horizontal primordial magnetic field, 
whereas galactic dynamo action may even be enhanced by 
buoyancy (Parker 1992; Moss, Shukurov & Sokoloff 1999); 
however we neglect the effects of magnetic buoyancy here 
so as to keep the model as simple as possible. 

Thus, any horizontal external magnetic field will be 
unimportant. However, the behaviour of the vertical mag- 
netic field B z is different. For a mean velocity field (in- 
cluding rotation) that is independent of the vertical coor- 
dinate z, a uniform B z is not distorted by the flow and 
is not subject to vertical diffusion (since d 2 B z /dz 2 = 0). 
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Therefore, such a field can be dragged into the disc by a 
radial flow. The dragging is especially efficient if the disc 
is surrounded by vacuum (where magnetic diffusivity is ef- 
fectively infinite) . The pile- up of the vertical magnetic field 
will be opposed and balanced by the radial magnetic diffu- 
sion. The resulting steady state was discussed by Lubow et 
al. (1994). A z-dependent rotation would add an axisym- 
metric azimuthal magnetic field produced from B z by the 
vertical shear. In addition, an azimuthal field can be pro- 
duced by differential rotation from a radial field generated 
by a vertical gradient in the inflow velocity. 

This picture is very much enriched by any dynamo 
action in the disc. The dynamo generates its own mag- 
netic field that can interact with the intruding vertical 
flux either constructively or destructively. The basic dy- 
namo mode in a galactic disc (where the a-effect is posi- 
tive in the northern hemisphere) is of quadrupolar parity 
(see, e.g., Parker 1979, Zeldovich et al. 1983, Ruzmaikin 
et al. 1988), opposite to the parity of the (quasi-) uniform 
vertical magnetic field. Therefore, the two fields can only 
interact via nonlinear effects. The situation can be differ- 
ent in accretion discs where magneto-rotational instability 
can produce a negative a-effect resulting in a predomi- 
nantly dipolar dynamo field (Brandenburg, Nordlund & 
Stein 1995). The synergy of the vertical flux and the dy- 
namo may result in physically nontrivial effects such as 
the semi-dynamo of Drobyshevski. Some of the behaviours 
obeserved in our models for large inflow speed can be con- 
sidered as a nonlinear analogue of semi-dynamo action 
wherein a subcritical dynamo action modifies non-trivially 
an invading external magnetic field. 

Radial flows can affect the disc dynamo in two obvious 
ways, both suppressive. As discussed by Moss, Shukurov & 
Sokoloff (2000), the radial advection hinders dynamo ac- 
tion, basically because the advection distorts the field from 
the eigenfunction. This effect can be important when the 
radial velocity u r <^ 2r)/h ~ lkms , for parameter val- 
ues typical of the Solar neighbourhood (h ~ 500 pc is the 
scale height of the ionized layer, rj ~ 10 26 cm 2 s~ 1 is the 
turbulent magnetic diffusivity). The second effect is due 
to the contribution of an imposed magnetic flux, dragged 
by the flow, to the nonlinear saturation of the dynamo ac- 
tion, customarily modelled by the a-quenching. Now the 
strength of the a-effect, responsible for the dynamo action, 
is reduced by the magnetic field as a oc (1 + B"* /Bq)~ x , 
where B is the strength of the local mean magnetic field, 
comprising both the dynamo and advected imposed flux. 
Here Bq is a reference field strength at which nonlinear 
dynamo effects become pronounced. The suppression of 
the dynamo action by the intruding magnetic field can be 
expected to be significant if it becomes comparable to Bq. 

All these effects cannot easily be separated, and it is 
difficult to assess their possible importance without the 
detailed analysis attempted here. In particular, our results 
presented below show that the dynamo can sometimes sur- 



vive a relatively strong imposed vertical flux - the dynamo 
can fight against its suppression. 

As discussed by, e.g., Lacey & Fall (1985) and Moss 
et al. (2000), the main mechanisms resulting in angular 
momentum transport with ensuing systematic radial flow 
in spiral galaxies are (in order of importance) the non- 
axisymmetric gravitational field of a bar or spiral arms 
(Lubow, Balbus & Cowie 1986), viscous and magnetic 
torques (Mestel 1999), and the infall of matter onto the 
disc (Pitts & Tayler 1996). The magnitude of the radial 
velocity and its radial profile are uncertain; the plausible 
values range from 1 to lOkms -1 at a distance of order 
3-10 kpc from the centre in a spiral galaxy such as the 
Milky Way. 

We examine two basic galaxy models. One has a gas 
density distribution taken from Milky Way observations in 
r > 2 kpc (see Sect. 3.2) but, in order to allow examina- 
tion of a wider range of parameters we limited the central 
gas density in r < 2 kpc to unrealistically small values. 
These calculations are described in Sects. 5.1 and |5.2| . In 
Sect. ^ we take a more realistic value of the central gas 
density, which is very strongly peaked in r < 2 kpc, and 
present results for a limited range of parameters applicable 
to the Milky Way. Thus readers who are only interested 
in applications to the Milky Way might consider omitting 
reading Sects. 5.1 and 5.2. 



We emphasize that we consider a prescribed dipolar 
flux, trapped by the disc at the time of galaxy formation. 
In reality, a galaxy is likely to trap a mixed parity flux, 
but the quadrupolar component will merely act as a seed 
for the dynamo. 

3. The model 

3.1. Basic equations 

We study a dynamo active disc embedded in a passive, 
conducting halo, and solve the standard a 2 to dynamo 
equation in the disc, 

QB 

— = V x (v x B + aB) - V x (tjV x B) , (1) 

but introduce an inwardly directed axisymmetric radial 
accretion flow and an axisymmetric vertical magnetic flux. 
The a-effect and turbulent magnetic diffusion are specified 
via scalar fields a and rj. We study axisymmetric solutions 
of Eq. (|l|) in the domain — Z < z < Z, r m j n < r < R, and 
define the aspect ratio of the domain as £ = Z/R;, where 
r and z are cylindrical polar coordinates. We choose 



ao sin(7rz//i) 



a = < 1 + B 2 /4irkpv 
, 



\z\<h, 
\z\ > h , 



defining the semi-thickness of the dynamo active disc, h, 
assumed to be independent of r (of course, h < Z). Here 
ao is the background (unquenched) magnitude of the a- 
effect, the gas density is denoted p (a function of r), and 
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v t is the turbulent velocity. We also define the disc aspect 
ratio A = h/R;. We have thus introduced a conventional 
a-quenching nonlinearity, with k a constant of order unity 
(reflecting some of the uncertainty in our ideas about the 
reaction of the large scale magnetic field on to the tur- 
bulent motions, but ignoring ongoing controversies about 
the nature and effectiveness of the feedback of the dynamo 
generated fields on the cv-effect - see, e.g., Vainshtein & 
Cattaneo 1992, Blackman & Field 2000, Kleeorin et al. 
2000, Brandenburg & Subramanian 2000). 

Since we will consider only axisymmetric solutions, we 
can write 



B = B^ + Vx , 



(2) 



where B^ and A^ are the azimuthal components of the 
magnetic field and vector potential, respectively. We also 
put 



= u(r)r + rtt(r)<j> + u dla 



(3) 



where Mdia = allows for the turbulent diamagnet- 

ism (Zeldovich 1956, Roberts & Soward 1975), u is the 
prescribed radial inflow, and Q is the galactic angular ve- 
locity. 

For the magnetic diffusivity, we set 



V 



i] , \z\< h 

770 [1 + (r?i - l)0(z)] , \z\>h 



where rjo is a constant, 

{\z\-hf 



O(z) = 1 — exp 



/-' 



is a smoothed step function and, rather arbitrarily, l n = 
— A). Thus, the constant rji is the ratio of the diffusion 
coefficient high in the halo to that in the disc. We consider 
models with the turbulent magnetic diffusivity in the halo 
significantly larger than in the disc. This appears to be a 
reasonable assumption for a galactic disc surrounded by 
a hot, turbulent halo (Sokoloff & Shukurov 1990; Poezd, 
Shukurov & Sokoloff 1993) as well as for an accretion disc 
corona supported by turbulent pressure (Ouyed & Pudritz 
1997). Anyway, the larger is the magnetic diffusivity out- 
side the disc, the easier is the dragging of the vertical 
magnetic flux. 



3. 2. A model for the Galaxy 

We take Q, from a model of the Milky Way (Burton & Gor- 
don 1978; see also p. 123 of Ruzmaikin ct al. 1988), and 
show V = rfl(r) in Fig. |l|. The surface density of the in- 
terstellar gas a is taken from Sanders, Solomon & Scoville 
(1984), and also shown in Fig. [l]for r > 2kpc. According 
to Sanders et al. (1984) (see also Dame 1992), the surface 
density of molecular hydrogen increases to 64 Mq pc~ 2 at 
r = 1 kpc and a few hundred Mq pc~ 2 inside that ra- 
dius. In our calculations we use both a central density of 
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Fig. 1. The radial dependence of the Galactic rotation 
velocity V — rfl (solid) (Burton & Gordon 1978), the ra- 
dial velocity in the disc midplane uo (multiplied by —50), 
Eq. (|4|) (broken), and gas surface density a at r > 2 kpc 
(Sanders et al. 1984) (histogram). 

300 M & pc~ 2 and also a naive extrapolation of Fig. ^ to 
5 Mq pc -2 at r — 0. The latter value is used to allow 
exploration of general properties of dynamos with inflow 
without the necessity of using a very high radial resolu- 
tion. These distributions of p and rQ describe a rather old- 
fashioned model of the Milky Way, with the Solar galac- 
tocentric distance and rotation speed being 10 kpc and 
250 km s , rather than the currently accepted 8.5 kpc and 
220 km s -1 , but we do not believe that this choice affects 
our conclusions. We neglected any possible z-dependence 
of f2 outside of the disc. This is rather uncertain, and 
would not dramatically affect the computed fields, given 
that the inflow and dynamo action are confined to the 
disc region, and so also is the bulk of the magnetic field. 
We neglect the possibility of any dynamo action in the 
halo (Sokoloff & Shukurov 1990; Brandenburg et al. 1992, 
1993). 

With 1] growing with z, turbulent magnetism, resulting 
in effective magnetic field transport at a velocity ttdia, 
compresses magnetic field to the disc (Brandenburg et al. 
1992; Gabov, Sokoloff & Shukurov 1996). However, this 
effect is relatively unimportant. 

Moss et al. (2000) discuss a crude but plausible model 
for the radial profile of a radial flow in the disc plane of a 
spiral galaxy, valid over all but the innermost part of the 
disc plane, 



Uoir ) = - 1A ] ^ e (r-Ro)/ro 

s 



M 



lM yr-i 



(4) 



where we put ro = 5 kpc. With the density profile shown 
in Fig. [I], the resulting mass accretion rate M ~ 2irrauo « 
lMQyr -1 is independent of radius within 20 per cent in 
5 < r < 15 kpc, and decreases at smaller r (mimicking a 
mass sink due to star formation in the galactic molecular 
ring). In order to avoid large flow velocities in the inner 
part of the disc, where observational data are poor and 
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the situation uncertain, we caused uq to go smoothly to 
zero between r — r u = 1.5 kpc and r = 0. Then if uq = 
— lkms _1 atthe Solar position, we have uq ks — 0.6kms _1 
at r = 300 pc, which corresponds to M w 0.3M©yr _1 at 
that radius (for a = 300M Q pc~ 2 ), as given by Morris & 
Serabyn (1996). We also explored the effect of reduction 
in the value of the truncation radius r u (see Sect. |6.2| ). 

The resulting radial velocity profile uq is shown in 
Fig. [I]. We also define a function 

1, \z\<h, 
exp[—(\z\/h - l) 2 ] , \z\ > h , 



fu(z) 



and put u(r,z) — uo(r)f u (z), ensuring that the inflow is 
essentially zero away from the disc. 

We consider a range of models with various values of 
Uo(Rq) in Eq. (Q) by replacing the factor 1.4 by 0.8U 
and varying Uq\ u (Rq) ~ — lkms -1 seems to be the 
most plausible value. 

3.3. Introduction of dimensionless variables 

Equation ([l]) is nondimensionalized in terms of a length 
R and a time h 2 /t]q. We measure angular velocity, radial 
velocity and a in units of their typical values Oo, Uq and 
cvo respectively, and B in units of B a = (inkpov 2 ) 1 / 2 , 
an equipartition value with respect to the turbulent ki- 
netic energy. The choice po = 5 x 10 _24 g cm -3 corre- 
sponds to the maximum surface density given by Sanders 
et al. (1984) outside the central region, which is about 
28M pc _2 , with the effective scale height of the neutral 
gas of 200 pc taken to be independent of r. The gas density 
p was taken to be independent of z in \z\ < h; its value in 
\z\ > h is irrelevant, as a — there, and p only appears 
in the a-quenching term. We took fio = 25kms -1 kpc~ 1 . 

The poloidal and toroidal parts of Eq. (Q) become re- 
spectively 
dA 

— = (\R a aB + \R u ur x B + A 2 w d ia X B) • d> 
at x ' 



\ 2 nV 2 A , 



(5) 



dB ^ 
— = V x (\R a aB + \R u ur x B + i? w ftr</> x B 
ot 

+ A 2 w dia x B) ■ 4> + X 2 vV 2 B , (6) 

where A and B are the dimensionless forms of A^ and B^ , 
and 

^ 2 9 2 1 d 1 d 2 

V 2 = 1 1 . (7) 

dr 2 r dr r 2 dz 2 

The radial velocity in Eqs. (|5|) and (^) is measured in units 
of Uo, giving 

u = -Omir/R)- 1 exp(r/r ) . (8) 

In these equations, and below, we refer to dimension- 
less quantities, unless otherwise stated. Thus we have the 
dimensionless parameters 



Ra 



no 



n h 2 



Ru = 



Uoh 



We chose R — 15 kpc, Z = 5 kpc, h = 500 pc (the scale 
height of the ionized gas, distinct from that of the neu- 
tral gas), i]i = 30, and considered two values of the 
disc magnetic diffusivity, 770 = 10 26 and 3 x 10 26 cm 2 s _1 . 
The former (more conventional) value of 770 determines 
R u = 18.75, the latter R^ — 6.25, and both suggest that 
if ao is of order 1 kms -1 , then values of R a of order unity 
are plausible; below we take R a = 0.5 and 2, respectively, 
both somewhat supercritical values. With 770 = 10 26 and 
3 x 10 26 cm 2 s _1 , values of the radial magnetic Reynolds 
number R u = 3 and 1, respectively, correspond to inward 
velocities of about u Q = — lkms -1 at the solar radius of 
10 kpc. With this value of uq(Rq), the unit radial speed 
is Uq ~ 1.8 kms" 1 . If vt = 10 kms -1 (at all radii), the 
unit of magnetic field is Bq = 8/c 1 / 2 p,G, and the unit of 
magnetic energy is ^BqR 2 Z. 

3-4- Boundary conditions 

At the boundaries z = ±£, we used the boundary condi- 
tions dB/dz — 0, B r = —dA/dz = 0. If the inner radius of 
the computational domain is at r m j n = 0, then regularity 
demands A = B = atr = r m ; n . For the many of our com- 
putations we took, for numerical convenience, r m j n = 0.02, 
and applied boundary conditions on A and B that were 
consistent with the analytical behaviour as r — > 0; most 
importantly dB z /dr = 0. We also used r m i n = for some 
computations - comparison between solutions with these 
values showed good agreement. The boundary condition 
at r = 1 is B = 0, and the initial field is a small amplitude 
seed of mixed parity. 

We set A — \B at all times on the outermost radial 
grid line r = 1 with B = const, and solve the difference 
equations resulting from Eqs. (||) and @ in the region 
Tmin < t < 1. With B z = r~ 1 d(rA)/dr, this specifies a 
constant vertical magnetic flux through the disc midplane 
z = 0, given by F = 2n J Q B z rdr = nB in dimensionless 
units. This prescription means that at time t = 0, when A 
is negligible at the inner grid lines, there is formally a dis- 
continuity in A, and so an infinity of B z , near the bound- 
ary r = 1. However this is immediately smoothed by dif- 
fusion, and the field is smooth after a few time steps. Our 
interest is mainly with stationary magnetic configurations 
(either steady or oscillatory) established at later stages 
of evolution, especially as initial conditions are physically 
uncertain. We confirmed, by performing a trial integration 
with the boundary condition that B z was uniform on the 
boundaries z — ±£, for a prescribed value of the flux, that 
our computational procedure led to the unique solution. 

Thus the interpretation of B is that it is the uniform 
field in the z-direction that would have the same total 
flux through the disc plane z = 0, r < 1 as the model 
considered. 

Equations (|5j) and (^|) are discretized on a grid with N r 
points uniformly distributed over the range r mm < r < 1 
and N z points distributed over — £ < z < £. The spatial 
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discretization is second order accurate in space, and a sec- 
ond order Runge-Kutta method was used to advance the 
solution in time. Our standard resolution was N r — 51, 
N z — 201, but we used N r — 401 in the models with large 
central density, cr(0) = 300 Mq pc~ 2 . 

3.5. Diagnostic quantities 

Four key parameters are used to quantify the properties 
of our stationary solutions, denoted by E, P, p and /. 
E is defined to be the total magnetic field energy in the 
computational domain. 

The global parity P is defined as P = (E cvcn — 
E dd)/E, where P C ven and P dd are the global energies of 
the parts of the field possessing respectively quadrupole- 
like and dipole-like symmetry with respect to z = (cf. 
Brandenburg et al. 1992). Thus a purely quadrupole-like 
field has P = +1, and P = — 1 for a dipole-like field. We 
similarly define Pi, the field parity in the region interior 
to f; = 0.02 (300 pc in dimensional units). We monitor P 
to allow for the possibility that dipolar magnetic field can 
dominate in the centre of the disc (with a predominantly 
quadrupolar global parity), as possibly observed in the 
Milky Way. As the dominant toroidal component of the 
magnetic field is strongly concentrated to the disc region, 
values of E, P and P calculated for the disc region will be 
close to the values calculated for the complete disc+halo 
system. 

The vertical magnetic flux within radius r is F(r) = 
2tt J q B z rdr — 2irA(r)r, whereas the total magnetic flux 
through the disc midplane is F = F(l) = irB. Then 
Pi = P(0.02) is the flux within radius r = 0.02 (aris- 
ing necessarily from the odd-parity part of the field) , and 
the ratio ip — F-JF is the fourth descriptor of our results. 
Note that for a strictly uniform field in the z-direction we 
have A(r) oc r, and ip = rf = 4 x 10~ 4 . 

To summarise, our results are presented in terms of 
the total magnetic energy E, the global magnetic parity 
P, the parity Pj of the field in the inner 300 pc of the 
disc, and tp, the relative vertical magnetic flux through the 
inner 300 pc. The average dimensionless vertical magnetic 
field within the inner 300 pc can be obtained as P 2jC = 
2.5 x l0 3 Bp. 

4. The reference states 

In this section we briefly describe two basic states of the 
system considered - firstly, dragging of an imposed verti- 
cal flux without any dynamo action and, secondly, a disc 
dynamo without any imposed flux or radial flow. These 
states provide references against which the relative im- 
portance of dynamo action and advection in the composite 
model can be assessed. For the dynamo model described 
in Sect. 4.2 we take the lower of the values discussed for 
the central gas density, a(0) = 5M Q pc -2 . The corre- 
sponding dynamo model for the high central density case, 




Fig. 2. The poloidal field lines (equally spaced in poloidal 
magnetic flux rA) in a stea dy s tate, for the kinematic 
advection calculation of Sect. |4.l| with R u = 18.75, ?/o = 
10 26 cm 2 s- 1 R a = 0, and (a) R u = 1.5, and (b) R u = 6. 



cr(0) = 3OOAf pc" 2 is described in Sect. EU|. The gas 



density does not enter the purely advective-diffusive cal- 
culation discussed in Sect. 



4.1 



4-.1. Magnetic flux captured by a dynamo inactive disc 

We start by considering the advection of imposed verti- 
cal flux by a radial inflow in a differentially rotating disc 
in the absence of dynamo action. Without the a-effect in 
Eq. (0), the equations (]5|) and ([|) for A and B decouple 
and their solutions as t — > oo can differ from zero only 
because 8^0. So it is clear that, for given values of R u 
and R u , solutions will scale with the single parameter £>, 
with magnetic energy P a( jv oc B 2 . Equilibrium is deter- 
mined solely by a balance of advection and diffusion, and 
the toroidal field B is parasitic on A, via twisting of the 
poloidal field lines by the differential rotation and vertical 
shear in u. The field parity is odd, P = —1, as this is the 
parity of the imposed external field. 

We show in Fig. || the poloidal field lines for tjq — 
10 26 cm 2 s~\ P w = 18.75, and R u = 1.5 and 6. The effects 
of increasing the magnetic Reynolds number of the radial 
flow are clearly seen: field advection is markedly stronger 
when R u is larger. 
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Fig. 3. The poloidal field lines (a) and contours of toroidal 
field (b), equally spaced in poloidal flux and toroidal field 
stren gth respectively, for the pure dynamo calculation of 
Sect. |J with R u = 0, R u = 18.75, r) = KFcn^s" 1 and 
R a = 0.5. The field has an even (quadrupolar, P = +1) 
structure with respect to the midplane (z — 0), so only 
the upper half space is shown. 



4-2. Disc dynamo without inflow 

When R a is non-zero, but R u — and B = 0, we have 
a standard disc dynamo. For reference, we show the field 
structure of a nonlinear solution with our standard pa- 
rameters R a = 0.5, R u = 18.75, 770 = 10 26 cm 2 s" 1 
in Fig. 0. The field parity is, as expected, quadrupolar 
(P = +1). The toroidal field is strongly concentrated in 
the disc, \z\ < 0.0333, where the ev-effect is non-zero, but 
the poloidal field pervades the halo. In this case we set 
B = but, in the absence of an inflow, the value of B is 
almost irrelevant. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
r 

Fig. 4. The poloidal field lines (a) and the contours of 
toroidal field (b), both equally spaced, for the calculation 
of Sect. |r| with R u = 3, R u = 18.75, ?y = 10 26 cm 2 s-\ 
R a — 0.5, B — 10~ 3 . Continuous and broken contours 
refer to positive and negative values respectively. 



fected by changes in a few key parameters. Results for our 
'standard' magnetic diffusivity are presented in Sect. 5.1 
and those using an enhanced diffusivity are in Sect. 5.2 



5. Dynamos with inflow; field advection and 
semi-dynamo action 

In this section we discuss axisymmetric steady state mag- 
netic fields in a dynamo active disc where trapped vertical 
magnetic flux is dragged in from large radii by an axially 
symmetric flow. We explore ranges of R u and B, the two 
less certain parameters of the model. The importance of 
the external magnetic field is also sensitive to the mag- 
netic diffusivity 77 in the disc - the dragging is more effi- 
cient when the magnetic field is more strongly coupled to 
the flow, i.e. when 77 is smaller for a given inflow speed. 
We consider two values of 77, one being a standard value 
for galactic discs, and the other a factor of three larger. 

We first consider the 'low' value for the central density, 
a = 5 Mq pc~ 2 , in order to allow the use of a lower radial 
resolution. In Sect. 5.1 results from a basic model, where 
we choose a single set of typical Galactic parameters, are 
given. Then we go on to explore how the results are af- 



5.1. Standard magnetic diffusivity: 770 = 10 2 

We begin by describing a typical nonlinear model with 
?7o = 10 26 cm 2 s , choosing parameters B = 10 -3 (cor- 
responding to a mean vertical field of 8 x lO^fc 1 / 2 G in 
imensional units), R u — 3, R a — 0.5 and R u — 18.75. 

We display in Fig. ^ the steady state magnetic field 
structure found in this case. There are clear deviations 
from pure parity, either odd or even. The global parity 
parameter is P — 0.85 indicating an almost quadrupolar 
parity, whereas in the inner 2% by radius P\ = 0.02. A 
small value of |Pi| means that the magnetic field is rela- 
tively small either above or below the disc, in this case in 
z > - see Fig. ||b. 

Our results are summarized in Table [l]. The fields 
evolve to steady configurations, except with the smaller 
B values with R u = 6, when the dynamo is barely excited 
when B — and oscillates for B = 0, 10 -3 and 10~ 4 . 
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Table 1. Summary of results for 770 = 10 26 cm 2 s _1 , the central gas surface density <r(0) = 5M Q pc~ 2 and R u — 
0, 1.5, 3, and 6. The two parity values given are respectively the global parity P and that in the inner 2% by radius, 
P. The values of E for B — are those for the pure dynamo, Pdyn- The magnetic energy obtained in the absence of 
dynamo action, Padv, scales with B 2 , but is given explicitly for convenience, ip = Fi/F is the fraction of the magnetic 
flux within radius r = 0.02, and B z ^ c is the dimensionless mean vertical field through z — within radius r — r\ = 0.02 
(corresponding to the inner 300 pc in radius). There are no entries for tp when B = 0, as there is zero flux through the 
midplane for a purely even parity field. Where a range of values are given, the solutions are oscillatory, and the values 
are maxima and minima. Note that the flux parameter B is the strength of the (dimensionless) uniform vertical field 
that would have the same flux through the midplane as the corresponding model. For reference, when R u = 0, B = 0, 
we find (P, Pi,E) = (+1, +1, 0.64). 
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R u = 3 
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4.5 x 10~ 2 




£>z,c 


52.5 


5.3 


0.55 
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4.10 
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0.125 
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0.021 


0.0067 . . . 0.030 


0.0051... 0.030 




£>z,c 


312.0 


29.3 


0.53 


0.017. . .0.075 


0.013... 0.0075 






5.1.1. Field parity and strength 

Two trends are apparent. As illustrated in Fig. |^, for given 
R u , increasing B causes the parity to become increasingly 
negative - i.e. more dipolc-likc. However, this effect is 
uniform neither in radius, nor between the poloidal and 
toroidal field components. We show in Fig. || the evolu- 
tion of the parity of the poloidal and toroidal parts of the 
field for the calculation illustrated in Fig. [|. Of course, as 
the energy in the toroidal part of the field is much greater 
than that in the poloidal part, the parity of the toroidal 
field dominates the global parity. For a given B and R u , the 
magnetic field in the inner disc r < 0.02 is more dipole-like 
than in the rest of the disc as a stationary field configura- 
tion is approached (i.e., Pi is systematically smaller than 
P). This also applies to the lower values attained by P 
and p in oscillatory solutions. 

Secondly, increasing R u decreases the dynamo effi- 
ciency, both in the linear and nonlinear regimes. This is 
masked to some extent when B 7^ by the increasingly ef- 
ficient inward advection of the vertical field. When B = 0, 



we find that the total dynamo-generated field energy de- 
creases as R u increases: the trend is consistent with the 
findings of Moss et al. (2000) who, however, mostly consid- 
ered linear solutions. We note that the magnetic energy 
(mean taken when oscillatory) of these two-dimensional 
solutions with B — varies approximately linearly with 
Ru as 

Pdyn «P (l-P«/6) , (9) 

with E = 0.64 for 77 = lO^cms -1 , < R u < 6 (see 
Table Q) and E = 0.90 for 77 = 3 x 10 26 cms" 1 , < R u < 
5 (see Table ||; we also note that Pdyn = 0.46 for R u = 3 
consistent with this fit). 

The result that advection inhibits dynamo action ap- 
pears quite general; Moss et al. (2000) produced an an- 
alytic argument, and our numerical results arguably ap- 
ply to an unfavourable the advection tends to 
sweep the field towards the central regions where the lo- 
cal growth rate is largest. 

The dynamo and the inward advection of the vertical 
field combine in a non-additive manner. In Fig. [?] (see also 
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Fig. 5. The field parity in the whole disc, P (solid), and 
in the inner region of the disc r < 0.02, P; (broken), ver- 
sus the vertical flux parameter B for the calculation of 
Sect. |5 1 1| with R u = 1.5, R u = 18.75, r\o = 10 26 cm 2 s-\ 
and R a = 0.5. Panels (a), (b), (c) are for the cases 
= 1.5,3,6 respectively. In (c), the range of variation 
of the oscillatory solutions is indicated when B < 10~ 2 . 
The field approaches purely odd parity, i.e. P — > — 1, as 
B increases, but the field near the centre of the disc has 
a systematically smaller value of P (i.e., is more dipole- 
like) than that in the disc as a whole, and needs a weaker 
external field B ~ 10~ 2 to become approximately dipolar. 




-0.5 



Fig. 6. The parity of the poloidal field (solid) and toroidal 
field (broken) in the whole computational domain as a 
function of time for the calculation illustrated in Fig. [|. 
(Time is measured in units of 7.5 x 10 8 yr). 
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Fig. 7. The total field energy E, versu s the flux param- 
eter B for the calculation of Sect. |5.1.l| . E is shown as a 
solid line for (a) R u = 1.5, (b) R u — 3 and (c) R u = 6. 
The other parameters are R u — 18.75, 770 = 10 26 cm 2 s~ x , 
and R a = 0.5. The solution for an advected field with no 
dynamo action, i.e. R a = 0, is shown by a broken line 
(cf. Sect. 4.1), and that for a dynamo without inflow, i.e. 
R u = 0, is shown by a dotted line (cf. Sect. fl[|). The field 
approaches that found in the absence of dynamo action 
as B increases, and that of a pure dynamo (without in- 
flow) as B decreases. Some solutions shown in panel (c) 
are oscillatory, and then the solid and dotted lines show 
the range of the field variation over a period. 



Table [j]), we show the magnetic energy in a stationary 
state, together with that resulting from pure advection 
(R a — 0), and pure dynamo action (B = 0). The total 
hybrid field energy E consistently exceeds _E a d v and is 
close to Edyn when B < 10 ~ 2 , whereas E < i? a dv for 
B = 1; we deduce that E — > -E a dv as B — > 00. 

The behaviour of the system when R u — 6 differs 
somewhat from that for the smaller values of R u . Now 
an oscillatory dynamo is excited when B = (see Fig. ||). 
The stationary state has relatively small energy and i? a dv 
exceeds Ea yn for B j£ 10 -3 . Nevertheless, dynamo action 
plays a significant role; in the range 10~ 3 < B ^ 1 dynamo 
action again reduces the magnetic energy of the system be- 
low that of the purely advected field, so that E < E^y. 
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Fig. 8. Upper panel: Magnetic energies in the odd (solid) 
and even (broken) parts of magnetic field in the whole 
computational domain for ijq — 10 26 cm 2 s _1 , B = 10 -3 
and R u = 6 (see also Table ^). The solution is oscilla- 
tory but its initial stage of growth is similar to that of 
non-oscillatory solutions. Lower panel: Global parity of 
magnetic field versus time. The even part of the magnetic 
field oscillates at a much larger amplitude than the odd 
part. Time is measured in units of 7.5 x 10 s yr. 



This case shows some of the features of the model with 
3 x 10 26 cm 2 



7/0 = a x iir" cm' s and R u — 10 discussed in Sect. 5.2. 

The oscillatory behaviour of the solutions at large R u 
appears in the linear solutions with B = also. However 
the approximate linear solution of Eqs. (||) and ||) devel- 
oped by Moss et al. (2000) in a one-dimensional approx- 
imation (retaining only z-derivatives, see their Eq. (Al)) 
remains non-oscillatory at any R u , and this is confirmed 
by their numerical solutions based on the 'no-z' approx- 
imation. Thus the oscillations appear to arise from the 
explicit fully two-dimensional structure of the model. The 
period of these oscillations is of the order of the galactic 
age, so they appear to be of little importance for galax- 
ies, but their origin in a broad dynamo theory context 
deserves further analysis. We have been unable to eluci- 
date further either analytically or physically the mecha- 
nism driving these oscillations. 

5.1.2. Field concentration towards the disc centre 

Another notable feature of the stationary states is that 
the fraction of the vertical magnetic flux stored in the 



central part of the disc, r < 0.02, exhibits a non-trivial 
dependence on both R u and B. Without dynamo action, 
ip, the fraction of the vertical magnetic flux stored within 
r = 0.02, increases with R u . These values of ip are identical 
to those shown in Table p] for B = 1, i.e. the solutions 
with large B are close to those for pure advection. Larger 
values of ip indicate stronger concentration of B z towards 
the disc axis. Note that <p = 4x 10~ 4 for a uniform vertical 
field. For dynamos without external vertical field (the last 
column in Table [j]), ip is not defined, because the field 
has quadrupolar parity and hence has no flux through the 
midplane z = 0. 

As illustrated in Table |l|, for R u = 1.5 (moderate in- 
flow), ip unexpectedly grows as B decreases below 0.1. This 
growth is less pronounced for a stronger inflow, R u = 3. 
The trend reverses as the inflow becomes very strong, with 
ip decreasing as B decreases. Thus, a larger trapped ver- 
tical flux results in a weaker central concentration of the 
vertical magnetic flux in the central part of a dynamo ac- 
tive disc provided the inflow in not too strong. 

For a fixed B, the field is more strongly concentrated 
to the centre when the inflow is stronger - a natural be- 
haviour. This behaviour occurs for B 10~ 2 at all values 
of R u explored and for R u < 6 for all values of B. Oscilla- 
tory solutions arising for R u = 6 behave differently having 
smaller maximum values of p than steady solutions with 
the same B but smaller i?„. 



5.1.3. Temporal variation 

Since the boundary conditions at r = R are inhomoge- 
neous when B ^ 0, Eqs. (JsJ) and @ do not admit expo- 
nentially growing solutions even at the linear stage when 
B 2 <C Ankpv 2 . However, the quadrupolar part of the mag- 
netic field, maintained by the dynamo, exhibits initial ex- 
ponential growth at an early stage followed by satura- 
tion, as might be expected for normal dynamo action. The 
dipolar dynamo modes decay at the values of R a and R^ 
considered here, so the dipole magnetic field within the 
disc only grows due to advection (perhaps modified by 
semi-dynamo action), and this growth is non-exponential. 
This is illustrated in Fig. || where we show the evolution 
with time of the energies in the even and odd parts of the 
magnetic field for R u = 6, B = 10~ 3 . The curve for the 
even field deviates from a straight line in this log-linear 
plot because of the contribution of the halo to the total 
energy (the exponential growth is confined to the disc). 
The solution is oscillatory, and the quadrupolar part has 
a much larger amplitude of oscillations than the dipolar; 
this is natural as the oscillations are definitely driven by 
the dynamo. The initial growth before the onset of the os- 
cillations is fairly similar to that found in non-oscillatory 
solutions. 

The time for the steady state configurations to become 
established will depend on the initial conditions. We have 
chosen what are arguably the least favourable, with no flux 
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initially penetrating the disc. The time to reach a steady 
state is then of the order of 10 10 yr. It is plausible that 
galaxies trap some poloidal flux at all radii as they form. 
Thus the timescale for our models to reach a steady state 
should be regarded as an upper limit. Moreover, changes 
in the disc, being directly driven by the inflow, will occur 
more rapidly if the inflow speed exceeds about lkms -1 , 
whereas the halo adjusts on a slower, diffusive, timescale; 
therefore, a stationary state can be established in the disc 
over a shorter time scale. 



5.2. Enhanced magnetic diffusivity: i)o = 3 x 10 26 cm 2 s 1 

We then studied a model with an increased magnetic 
diffusivity in the disc, ?yo = 3 x 10 26 cm 2 s — 1 and thus 
also in the halo. As in the previous section, we adopted 
a = 5 Mq pc~ 2 in r < 2 kpc. Now field amplification by 
differential rotation is weaker because of enhanced diffu- 
sion, R w — 6.25, and so the critical value of R a is sig- 
nificantly increased; we took R a — 2, again a somewhat 
supercritical value for dynamo action. With these param- 
eter values, R u = 1 corresponds to an inflow speed of 
about lkms -1 near the Sun. Our results are presented 
in Table |[ Comparing these with those in Table [j], we 
see a broadly comparable pattern in the results when 
Ru < 5, with advection significantly affecting the dynamo 
field when B <; 10~ 3 . Arguably, in this case advection 
is somewhat less effective, consistent with what can be 
expected for the larger magnetic diffusivity. In particular, 
the concentration of the vertical magnetic flux towards the 
disc axis, ip, increases as B decreases for 10 -4 < B < 1CP 1 , 
and decreases for smaller values of B. 

For R u = 10, the situation is rather different. Now 
the dynamo is not excited when B — 0. Even so, dynamo 
action plays a role; in the range 10 < B cxt < 1, we find 
that the saturated magnetic energy is less than that of the 
advected field, i.e. E < -Eadvj by a factor of approximately 
3 for the smaller values of B. The field always has a dipolar 
parity, P = — 1, not typical of the dynamo configurations 
with moderate dynamo numbers explored here. 

When a quadrupolar field is not excited, the dipolar 
field energy in the presence of an imposed flux and dynamo 
action is reduced from that when the dynamo is turned off. 
This leads us to deduce that a form of semi-dynamo, ex- 
citing an odd parity field, is operating (cf. Drobyshevski 
1977); i.e. subcritical dynamo action that is dependent 
on an external source of flux. The fields so generated, al- 
though of the same, odd, parity as the advected field, nev- 
ertheless interact destructively with the advected field, re- 
ducing the total energy below that of the purely advected 
field (i.e. that found with R a = 0). 

6. A Milky Way model with high central density 

Here we examine a model more specifically relevant to 
the Milky Way, with <r(0) = 3OOAf pc -2 and r] = 




Fig. 9. (a) The poloidal field lines equally spaced in 
poloidal flux and (b) contours of toroidal field with con- 
tour levels starting at unity and increasing by a factor of 
1.5 in toroidal field strength between contours, for the pure 
dynamo calculation of Sect. 3.1 with R u = 0, R^ = 18.75, 
i] = 10 26 cm 2 s- 1 , R a = 0.5 and cr(0) = 300M Q pc- 2 . 
The field has an even (quadrupolar, P = +1) structure 
with respect to the midplane (z = 0), so only the upper 
half space is shown. 



10 26 cm 2 s -1 . We take R u = 3, so the inflow velocity near 
the Sun is about 1 kms -1 . The reference model with no 

in this lin- 



dynamo action is as described in Sect. 4.1 



ear calculation the gas density has no role. We describe 



in Sect. 6.1 a nonlinear dynamo model with no radial 



flow, before giving our results for models with inflow and 
a range of B values. 

6.1. The basic dynamo model 

In Fig. ^| we show the poloidal field lines and toroidal field 
contours for the basic dynamo model (R u = 0, B = 0). 
Comparison with Fig. || shows that the toroidal field is 
much more strongly concentrated to the galactic centre 
than in the models of Sect. ||even though its magnitude at 
r ~ Rq remains approximately the same. This is a natural 
consequence of the strongly enhanced central density, and 
correspondingly increased field strength Bq at which the 
dynamo saturates. 

6.2. Dynamos with imposed flux through the midplane 

We first discuss results with R a = 0.5, the value used in 
Sect. 5.1 . In Fig. [l(] we present the field structure (poloidal 
field lines and contours of toroidal field) for a model with 
B = 10 -3 . Marked departures from pure symmetry are 
visible; here P = 0.87 and Pi = 0.82, although these 
are not as strong as for the corresponding model with 
small central density, see Fig. |]. The toroidal field is now 
much more strongly concentrated to the galactic centre 
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Table 2. Summary of results for calculations with 770 = 3 x 10 26 cm 2 s \ the central gas surface density cx(0) = 
5M pc~ 2 , and R u = 0, 2, 5, and 10. The format is similar to that of Table [|. When R u = 0, B = 0, we obtain 
{P,P l7 E) = (+1, +1,0.90). 
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160. 


4.50 


1.93 




0.21 




0.021 













Ru = 


10 












P 


-1 


-1 


-1 




-1 




-1 






Pi 


-1 


-1 


-1 




-1 




-1 






E 


7.0 x 10 3 


3.2 x 10 1 


3.0 x 10" 


-2 


2.6 x 10" 


-4 


2.6 x 10" 


6 


decay 


E a dv 


7.2 x 10 3 


7.2 x 10 1 


7.2 x 10" 


-2 


7.2 x 10" 


-4 


7.2 x 10" 


6 





<£_ 


0.19 


0.13 


0.032 




0.022 




0.022 






B z ,c 


475. 


32.5 


0.80 




0.055 




0.0018 








Table 3. Summary of results for 770 = 10 26 cm 2 s 1 , R u = 3, and the central gas surface density cr(0) = 300 Mq pc 2 . 
The format is similar to that of Table 0. For R u = 0, B = we obtain (P, P t ,E) = (+1, +1, 0.84) when R a = 0.5. and 
(+1, +1,0.17) for R a = 0.25. 



B 


1 


nr 1 


6 x 10~ 2 


3 x 10" 2 


2 x 10 -2 


10~ 2 


10~ 3 


10- 4 













Ra 


= 0.5 










P 


-1 


-1 


-1 


-0.86 


-0.69 


-0.35 


0.87 


0.990 


+1 


Pi 


-1 


-1 


-1 


-0.61 


-0.12 


0.07 


0.82 


0.998 


+1 


E 


7.66 x 10 3 


48.0 


4.6 


0.80 


0.88 


0.80 


0.44 


0.42 


0.42 


-E a( Jv 


7.68 x 10 3 


7.68 x 10 1 


27.6 


6.91 


3.07 


0.768 


7.68 x 10~ 3 


7.68 x 10" 5 







0.026 


0.021 


0.0089 


0.0032 


0.0051 


0.0080 


0.054 


0.067 




Bz,c 


65.0 


5.2 


1.3 


0.24 


0.26 


0.20 


0.14 


0.017 













Ra 


= 0.25 










p 


-1 






-1.0 


-0.996 


-0.56 


0.91 


0.999 


+1 


p 


-1 






-1.0 


-0.992 


0.12 


0.94 


0.999 


+ 1 


E 


7.64 x 10 3 






1.03 


0.38 


0.18 


0.11 


0.11 


0.11 


Eadv 


7.68 x 10 3 






6.91 


3.07 


0.768 


7.68 x 10~ 3 


7.68 x 10" 5 





£_ 


0.026 






0.0065 


0.0030 


0.0069 


0.025 


0.027 




Bz, c 


65.0 






0.49 


0.15 


0.17 


0.063 


0.0067 






with B\ r= o. 2/B\ r=Re ~ 20. Our results for < B < 1 
are given in Table ||. We s ee a broad similarity with the 
results of Sects. [O] and 5^2. The steady-state magnetic en- 
ergy decreases with increasing inflow rate, cf. Eq. (||). As B 
increases at fixed R u , then E^ ya < E < -E a dv and both P 
and Pi decrease monotonically, so the field becomes more 
dipole-like both globally and in the inner 300 pc. The only 



significant difference to the results found above with the 
smaller value of the central density is that Pi < P for in- 
termediate values of B, i.e. the central field is less dipole- 
like than the global field. Further (but related) we note 
that for these intermediate values of B, tp is smaller than 
for the corresponding model with cr(0) = 5 Mq pc -2 (see 
Table |l|) . We suspect this is because the stronger dynamo 
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netic configurations. This may be considered as another 



manifestation of semi-dynamo action, cf. Sect. 5.2 



-0.3 



Fig. 10. (a) The poloidal field lines equally spaced in 
poloidal flux and (b) contours of toroidal field with con- 
tour levels starting at unity and increasing by a factor of 
2.6 in toroidal field strength between contours, for the cal- 
culation of Sect. |6^ with R a = 0.5, R u = 3, = 18.75, 
r/o = 10 26 cm 2 s-\i? Q = 0.5, <r(0) = 300Af Q pc" 2 and 
B = 10~ 3 . Continuous and broken contours refer to posi- 
tive and negative values respectively. 

generated field now present in the central regions is better 
able to interact destructively with the advected field. 

We also experimented with a reduction of the trunca- 
tion radius r u for the velocity profile to 0.05 (correspond- 
ing to 0.75 kpc) from our standard value r u = 0.10 (see 
Sect. 3.2). This resulted in an increase of ip by a factor of 
about 2 when B = 10~ 3 , and of about 2.5 when B = 0.1. 
We consider this to give an estimate of the uncertainties 
in our model. 

Given the importance of the destuctive interaction of 
dynamo and advected fields, we also investigated a model 
with weaker dynamo action, taking R a — 0.25. These re- 
sults are also presented in Table ||. The change of the diag- 
nostic parameters as B increases follows a similar pattern 
as described above when R a — 0.5, except that the transi- 
tion to a dipolc-like field structure occurs at smaller values 
of B. In particular, we note that except when B = 3 x 10 -2 



(when already P ~ —1), the mean central field B zc is 
smaller than when R a = 0.5. We conclude that the dipolar 
magnetic field at the disc centre is supported by dynamo 
action even though the latter favours quadrupolar mag- 



7. Discussion and conclusions 

7.1. Implications for magnetic field in the Galactic centre 

An unexpected result of our simulations is that the rela- 
tive concentration of vertical flux in the central region, ip, 
does not increase monotonically with the vertical flux B. 
This effect is especially strong for moderate inflow speeds. 
For example, in Table || with R u = 3, R a =0.5, the frac- 
tion of the vertical flux stored in the inner 300 pc grows 
by a factor of 20 as the magnitude of the flux parameter 
B decreases from 3 x 10~ 2 to 10~ 4 . Increasing the cen- 
tral density enhances the central quadrupolar field, but 
reduces the dipolar vertical field there. 

Our model, with a high turbulent magnetic diffusion 
in the halo, favours advection (also the field is not an- 
chored high above the disc). Nevertheless, our results in- 
dicate that the nonlinear interaction with the dynamo lim- 
its the vertical magnetic field at the centre to relatively 
low values even when both the advection and the cap- 
tured vertical field are rather strong. In Tables [[], || and 
H we give the values of B ZjC = Btp/rf = 2500Bip, the di- 
mensionless mean vertical field within the inner 300 pc 
radius (i.e. r < n), the overbar denoting a horizontal av- 
erage. For the high central density model of Sect. ^ (Ta- 
ble |), the mean central field B zc is approximately con- 
stant in an interval 10~ 3 Ss B Si 3 x 10~ 2 . Quite generally, 
~B~z, c ^ (10 - 20)fc 1/2 ^G for B < 0.1 and realistic values 
of Ru, and generally decreases with B. (Note that this is 
a higher estimate than might be deduced from Fig. O, 
where B = 10~ 3 .) 

The estimated magnitude of B z ^ c may be further in- 
creased by a factor of 2-3 because of uncertainty in the 
inflow speed in the central region. However, any increased 
turbulent velocity in r < 300 pc (not allowed for in our 
models) would result in larger values of Bo, thus both en- 
hancing the quadrupole magnetic field due to the dynamo 
and increasing the turbulent diffusivity. The latter would 
hamper advection of the vertical magnetic field in the cen- 
tral regions, and could possibly compensate any increase 
in B z . c . Another possibility is that ao may be underesti- 
mated at small radii in our model if ao oc f2; this would 
make the dynamo stronger in the central regions, again 
favouring the generation of quadrupolar magnetic fields 
there. The resulting average vertical field at r < 300 pc 
would then remain not stronger than 0.01 mG, two orders 
of magnitude weaker than required if the magnetic field 
illuminated in the filaments fills the volume. 

We show in Fig. [ll] the radial variation of the az- 
imuthal and vertical magnetic field components in a model 
with a radial velocity of about lkms" 1 near the Sun, 
where the degree of concentration of B z to the centre, 
ip, is nearly the strongest among the models studied. 
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Fig. 11. Radial profiles of the azimuthal (a) and vertical 
(b) magnetic field components at z = — 500 pc (solid), 
z = (dashed) and z = 500 pc (dotted) in the inner 2 kpc 
in a model with R u — 3, R a — 0.5, rj — 10 26 cm 2 s , the 
central gas surface density 3OOM0pc -2 and B — 10~ 3 in 
dimcnsionless units. Magnetic field is in units of 8k 1 / 2 fiG, 
and radius in units of 15 kpc. 



The profiles are shown for z — and ±500 pc to pro- 
vide further indication (in addition to the parity param- 
eter Pi introduced above) of the degree of symmetry of 
the vertical profiles with respect to the midplane. Even 
though the dipolar component [with B$(z) — —B^—z) 
and B z {z) = B z (—z)] is significant (especially in the az- 
imuthal field), the dynamo- favoured quadrupolar compo- 
nent [with B<p(z) = B^(—z) and B z (z) = — B z (— z)\ is 
stronger. Correspondingly, the parity in the inner 300 pc 
is rather close to unity (see Table ||) . Further details for 
these high central density models are presented in Table [| 
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and in Fig. 

The vertical magnetic field shown in Fig. [Il] has a max- 
imum at r = 0; this is a typical feature. However, the 
maximum value is B z ~ 2 x 10~ 6 G at r, z = (taking 
k — 1), which is significantly below the level corresponding 
to equipartition with kinetic energy, B ~ 10~ 5 -10~ 4 G. 
The azimuthal field vanishes at r = (because of the az- 
imuthal symmetry) and reaches a maximum at a distance 
of about 500 pc from the centre; the maximum value of 
about 1.5 x 10 -4 G (again taking k = 1) is slightly above 
the equipartition level. 

To limit the range of B of physical interest, we can 
apply three constraints to our models. On observational 
grounds we can probably rule out solutions with E 



at 
o 



o 



OA 
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k,': 
o 




0.04 
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Fig. 12. (a) The total magnetic energy E versus the cap- 
tured flux parameter B in the model of Sect. ^ is shown 
by a solid line for R a — 0.5. The advected field energy 
and dynamo energy are shown dashed and dotted as in 
Fig. 0. (b) The mean vertical magnetic field (solid) and 
the relative flux ip through the inner 300 pc radius (bro- 
ken) versus B for R a = 0.5. The corresponding results for 
weaker dynamo action, with R a — 0.25, are shown with 
open circles and triangles, respectively. Magnetic field is 
measured in units of 8k 1 / 2 fiG. 



-Edyn, simply because Ed yn is close to the kinetic energy of 
the turbulence and corresponds well to the field strengths 
observed in the outer parts of the Milky Way and other 
galaxies. Thus larger values of the energy would imply 
strongly super-equipartition large scale (i.e. regular) fields 
throughout the galactic disc, for which there is no evidence 
in normal spiral galaxies (e.g., Beck 2000). This constraint 
thus suggests that we should probably restrict our atten- 
tion to values B 5s 5 x 10~ 2 . 

Secondly, there is evidence to suggest that, in common 
with external spiral galaxies, the global magnetic field of 
the Milky Way should be of even (quadrupolar) parity 
(Frick et al. 2000 and references therein). If conservatively 
we require P > 0, then we need B ^ 10~ 2 . 

Thirdly, if we estimate the maximum large-scale mag- 
netic field strength in the protogalactic gas to be ^ 3 x 
10 -12 G (Sciama 1994), then in contracting from 100 kpc, 
an estimated radius for the protogalaxy, to the present- 
day cylindrical radius of 20 kpc without loss of flux, we 
obtain a mean vertical field through the disc of less than 
about 10~ 10 G, i.e. B ^ 10~ 5 . This is a negligibly weak 
field in the present context. For a galaxy formed in a rich 
cluster of galaxies (though mostly populated by ellipticals 
rather than spirals), the external field could be a fraction 
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of 10 -6 G, corresponding to B ~ 0.1-0.01 in dimensionless 
units. This is clearly a very generous upper limit for the 
Milky Way. 

These considerations all suggest that we need only con- 
sider models from Table| with B <, 0.05 (of order 10" 7 G), 
and possibly only those with B Ss 0.01, or even much less. 
This implies that B ZtC 1.0, that is the mean vertical 
field near the galactic centre is less than about 10/c 1//2 /xG. 
Uncertainties associated with the choice of the inflow ve- 
locity cut off parameter r u (see Sect. 3.2) seem unlikely 
to change this estimate by more than a factor of 2 or 3. 
We also note that the same arguments applied to the low 
central density models of Table [j] with R u = 3 lead to 
similar estimates. 

Dynamo action is not the only obstacle to the accu- 
mulation of dipolar magnetic flux in the Galactic centre. 
Even the pure field advection models with the lower value 
10 26 cm 2 



?/o = llT"cm"s and R a = as in Sect. 4.1 (ignoring 
for a moment that these models have strictly odd parity, 
P = — 1 and so do not reproduce magnetic field structure 
near the Sun), and B 10~ 2 again yield B ZyC < 1, and so 
the mean central field is SS 10/iG. 

The models explored here indicate that a vertical mag- 
netic field of an average milligauss strength cannot be ac- 
cumulated in the centre of the Galactic disc. The striking 
dipolar symmetry implied by the straight shape of the 
magnetic filaments in the Galactic centre also does not 
develop in our models. Our model is not designed to pro- 
duce such fine filamentary structure. However, unless the 
filling factor of the magnetic filaments is extremely small 
(i.e. the filaments are not parts of a more uniform field illu- 
minated only intermittently - see Morris & Serabyn 1996 
for a discussion), there is certainly insufficient vertical flux 
present in our models. 

We emphasize that we have only considered magnetic 
fields resulting from Galactic-scale accretion and dynamo 
action. It cannot be excluded that the central regions of 
the Galaxy represent an autonomous system relatively 
weakly connected with the outer parts, at least in terms 
of magnetic field. In that case, it would be necessary to 
develop a separate, specific MHD model (allowing for dy- 
namo action) for this region. What we have shown, how- 
ever, is that the observed milligauss vertical magnetic 
fields within the inner 300 pc would be difficult to explain 
by the accumulation of magnetic flux from the outer parts 
if dynamo is active in the Galaxy, even if the total vertical 
magnetic flux is conserved. This conclusion seems rather 
general and independent of the actual dynamo mechanism. 

Our conclusions are in variance with those of Chan- 
dran et al. (2000), principally because these authors as- 
sume the large-scale magnetic field to be frozen into the 
interstellar gas. They argue that the dissipation of mag- 
netic energy in the hot gas, expected to occur at the proton 
gyroradius scale, can take as long as 10 14 yr. However, the 
large-scale magnetic field is more plausibly anchored in 
the warm phase (Beck et al. 1996) which has larger filling 



factor and can be pervasive, whereas the hot gas is buoy- 
ant and can leave the disc in a short time to be confined to 
isolated regions in the disc. Furthermore, the large-scale 
magnetic field is affected mostly by turbulent magnetic 
diffusion which is determined by the energy-range scale of 
turbulence about 100 pc, and the velocity at this scale, 
rather than by the dissipation scale of the turbulence. 
Moreover, mean-field dynamo action does not require gen- 
uine turbulence at all (i.e. any details of turbulent energy 
cascade are unimportant) - what is needed is merely ran- 
dom, mirror-asymmetric, flow. Just as fundamentally, our 
results indicate that dynamo action in itself acts to in- 
hibit advection of field: this is the reverse side of inflow 
inhibiting dynamo action in the absence of any external 
magnetic field. There is no reason to believe that this re- 
sult is dependent on the detailed nature of the dynamo 
(e.g. the a-effect). 

Chandran et al. (2000) similarly argue that gravity and 
tension of the large-scale magnetic field will suppress the 
random nature of interstellar gas motions thereby sup- 
pressing turbulent diffusion. It is, however, known that 
turbulent diffusion remains effective even in stronger grav- 
ity fields, e.g., in the solar convective zone. Magnetic fields 
undoubtedly suppress turbulent diffusion together with 
the a-effect (e.g., Krause & Radler 1980), but these effects 
plausibly occur simultaneously, so that dynamo action, 
controlled by the dynamo number D oc a/n 2 , can continue 
(e.g. Moss in Mestel 1999, p. 222; Brandenburg 2000.) 
Moreover, a large-scale magnetic field can even enhance 
turbulent magnetic diffusivity in some models (Branden- 
burg & Subramanian 2000). Efficient magnetic turbulent 
diffusion does not imply that gas will be removed from the 
disc together with magnetic fields (Chandran et al. 2000) 
as gas freely slips down the vertical sections of tangled 
magnetic lines. 

Altogether, we conclude that there are no compelling 
reasons to consider the large-scale galactic magnetic field 
to be frozen into the interstellar gas and the dynamo ac- 
tion to be suppressed. We reiterate that intense turbu- 
lent motions are observed in the interstellar medium, and 
these motions cannot avoid tangling the large-scale mag- 
netic field. 

7.2. More general implications for galactic dynamo 
theory 

Our system also demonstrates several features of general 
interest in dynamo theory. In the main they have been 
described above, and are merely recapitulated here. 

For large enough inflow velocities, the dynamo is killed 
in the absence of an imposed vertical flux, supporting the 
conclusions of Moss et al. (2000). As this situation is ap- 
proached with increasing R u , the dynamo may become 
oscillatory at relatively low energy; these oscillations ap- 
pear only in the fully 2D model studied here, and not in 
the approximate dynamo solutions discussed by Moss et 
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al. (2000). For small and moderate R u , the steady-state 
magnetic field scales as (1 — cRu) 1 / 2 with c close to 1/6. 

In the presence of a vertical flux (B^ 0), semi-dynamo 
action may occur; this is explicit when R u is large enough 
that pure dynamo action is suppressed, and the steady- 
state magnetic energy still differs significantly from that 
obtained with pure advection, i.e. without any cv-effect. 

In general the energy of the hybrid field (from dynamo 
action and advection) is not the sum of the energies of the 
corresponding pure dynamo and purely advected fields - 
for larger B values, E < E a( i v even in cases where a dy- 
namo is not excited when B = 0. This is somewhat remi- 
niscent of a result found by Moss (1996) when investigat- 
ing stellar dynamos with an imposed background field. We 
find that, as long as -E a dv <i Ed yn , the (true) dynamo gen- 
erated fields interact destructively with the advected field, 
i.e. resulting in E < E^dv In contrast, the interaction is 
constructive when i? a( jv < -Edyn, i-e. E > i? a dv + -Edyn in 
the stationary state. 

Although it seems plausible that a large enough verti- 
cal flux and inflow velocity could seriously inhibit dynamo 
action via cv-quenching caused by the advected vertical 
field, we do not see unambiguous evidence for this in our 
models. For large enough R u , the dynamo is in any case 
suppressed by essentially linear effects (even when B = 0). 
However, as B increases, there is evidence (E < -E a dv) that 
the advected field promotes some form of dynamo action, 
even in the extreme case R u = 10, B = 1 (Table ||). It is 
clear that the effects of nonlinearity in this dynamo system 
can take subtle forms. 

We have here discussed only dynamo systems with 
R a > 0, in which the basic field generated is of even 
(quadrupolar) parity, as believed to be appropriate to spi- 
ral galaxies. If in accretion discs R a is negative (Branden- 
burg et al. 1995, Riidiger & Pipin 2000), then a dipolar 
(odd) parity dynamo field is generated. The outcome of 
the interaction of such a dynamo with a trapped verti- 
cal flux of the same odd parity is unclear. We intend to 
address this question elsewhere. 
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